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Abstract 

In the present work it is studied the fermionic van Hemmen model for the spin glass (SG) with a transverse 
magnetic field T. In this model, the spin operators are written as a bilinear combination of fermionic operators, which 
allows the analysis of the interplay between charge and spin fluctuations in the presence of a quantum spin flipping 
mechanism given by F. The problem is expressed in the fermionic path integral formalism. As results, magnetic phase 
diagrams of temperature versus the ferromagnetic interaction are obtained for several values of chemical potential /i 
and T. The T field suppresses the magnetic orders. The increase of yu alters the average occupation per site that affects 
the magnetic phases. For instance, the SG and the mixed SG+ferromagnetic phases are also suppressed by fi. In 
addition, fi can change the nature of the phase boundaries introducing a first order transition. 
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It is well established that disorder can be a source of non-trivial effects in condensed matter physics. From 
Anderson localization to spin glass, there are plenty of examples supporting such role. As one more recent example 
there are indications that the presence of disorder is responsible by several complex effects in strongly correlated 
systems as, for instance, manganites or Cerium compounds (see Ref. (JJJ and references therein). The important point 
is that disorder can affect the charge sector as well as the spin one in such way that it becomes a constant cause of new 
interesting problems. From the theoretical point of view, the task of facing this kind of problems is surely a source of 
new approaches. 

Recent experimental findings in two specific Cerium compounds can illustrate the last comment in the previous 

CN ' paragraph. CePd\- x Rh x 12fl and CeNh- x Cu x 13[] are physical systems in which disorder combined with RKKY and 

Kondo interaction is responsible by the onset of a complex scenario. Their phase diagrams display a glassy state, 
ferromagnetism and a region dominated by the Kondo effect. For CeNi\~ x Cu x , several works have used the so called 
Kondo-Ising lattice (KIL) model |4J| to reproduce the experimental phase diagram of that compound. The Ising part 
of the model represents an intersite interaction between localized f-spin operators. It is precisely in the spin-spin Ising 
coupling Jjj that various types of randomness have been tested trying to represent the effects of the disorder in the 

f^i ■ CeNi\- x Cu x |5L la [i ia |9[] . For instance, the random coupling /y has been chosen as the Sherrington-Kirkpatrick 



(SK) model 1 1 QTI or as the Hopfield model fll ill . However, none of the choices previously mentioned were entirely 
satisfactory to describe how the glassy state in the CeNi\- x Cu x is replaced by an inhomogeneous ferromagnetism 
when the temperature is decreased yfl . 

Very recently, a new approach has been proposed for the problem discussed above, this approach uses the random 
coupling Jij as that one introduced in the van Hemmen (vH) model 11211 . This model was originally conceived to study 
the spin glass (SG) problem with classical Ising spin variables. The results indicated that the vH type of randomness 
in the KIL model is more adequate than those ones used in SK or Hopfield models to reproduce some aspects of the 



experimental phase diagram of the CeNi\- x Cu x 11311 . For instance, the previously mentioned ferromagnetic phase 
appears below the spin glass one. 

The improvement obtained using the coupling /y given as the vH model to study the CeNi\- x Cu x suggests the 
usefulness of this kind of randomness could be more general. In other words, it could also be used when it is present in 
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other quantum process than the Kondo effect. In this case, models with one type of fermions could be built by adding 
to the vH Ising term, for instance, a hopping, a pairing coupling 11411 or a transverse magnetic field. Nevertheless, 
surprisingly little consideration has been given to a fermionic version of the vH model. In fact, even less consideration 
is given to the fermionic vH model when a term capable to introduce quantum dynamics in the problem is added to 
this model. In fact, it lacks an investigation of the interplay of charge and spin fluctuations within a fermionic version 
of the vH model even for a simple quantum process provided by a transverse magnetic field. 

The purpose of the present work is to study the vH model in the presence of a transverse field F in a fermionic 
formulation which means that the Ising spin operators are written as bilinear combinations of creation and destruction 
fermionic operators. For this particular representation of the spin operators, the natural tool to obtain the thermody- 
namics is the functional integral formalism with Grassmann fields. A peculiarity of this representation is that S z t acts 
in a space with four eigenstates. Among them, there are two which are nonmagnetic 11511 . Therefore, we calculate the 
grand partition function for two situations. In the first one, here called AS fermionic van Hemmen (45-FvH) model, 
the four states are preserved. For this particular model, the chemical potential fi appears as an important parameter in 
the problem. In the second one, the 25 fermionic van Hemmen (25 -FvH) model, introduces a restriction, what allows 
considering only the two magnetic states IU5 U16I1 . One of the most important advantages of the vH model in both 



versions 4S and 2S is that the disorder can be treated without the use of the replica method. Therefore, the subsequent 
difficulties associated with this technique 11711 can be avoided in the present work (see discussion below). However, 



the problem is still treated within the static approximation (SA) in which the time fluctuations are neglected 1 18]. 

Indeed, the fermionic representation of Ising spin operators for the study of the SG problem is well established. 
For the SK type of randomness, the so called fermionic Ising spin glass (FISG) model Ilo.ll9ll20ll shows an interesting 
characteristic. Charge and spin correlation functions are connected. The average occupation of fermions per site and 
the replica diagonal element of the SG order parameter matrix are directly related. This relationship implies that, for 
instance, the nature of the PM/SG phase boundary can be modified by variations of chemical potential with the onset 



of a tri critical point [21[j. The presence of an additional transverse magnetic field F in the FISG model also introduces 



important changes in the PM/SG phase boundary. The first order part of that phase boundary is diminished when F 



increases [22]. Consequently, the position of the tricritical point is also affected. Therefore, one could expect that 
in the 4S-FvH model, variations of fi could also bring important consequences on the location and nature of phase 
boundaries. 

It should be remarked that the location of a first order PM/SG phase boundary is not a trivial task for models in 
which the replica method is used as the FISG one 122L 12311 . This difficulty has already been a subject of controversy 



for the classical Ghatak-Sherrington model (see, for instance, Refs. 1241 12510 . It is clear that this kind of controversy 
would not be needed in the FvH model. On the other hand, vH model (classical or fermionic) has a shortcoming. 
It lacks to present a multiplicity of metastable-states 12611 . Nevertheless, the vH model produces some results which 



are close to those ones found experimentally in SG physical systems (see Ref. 12711 ). Thus, the proper understanding 
of the vH model in a fermionic formulation can be important not only as a SG problem itself but also as something 
useful for applications as pointed out above and very recently discussed in Ref. 11311 . 

It should also be mentioned that the vH model with an additional transverse magnetic field F has already been 



studied. However, in a formulation where the classical Ising spin variables were replaced by Pauli matrices [28]. 
In this formulation one can observe important consequences for the phase diagram produced by the presence of T. 
The obtained phase diagram shows that the increase of F destroys the SG, the ferromagnetic (FM) and the mixed 
(FM+SG) phases found in the classical phase diagram. Nevertheless, this kind of formulation is not able to capture 
the consequences for the phase diagram when both charge and spin fluctuations are affected by F. 

This paper has the following structure. In the second section, the model is introduced and the grand canonical 
potential is derived. In the third section, several phase diagrams are constructed with particular attention for the 
effects of F and fi by solving the equations for the order parameters. The last section is dedicated to the conclusions. 

1. Model 

The Hamiltonian of the fermionic van Hemmem (FvH) model with a magnetic transverse field F is given by 
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s * = 2 [c W ~ i Ci ^' § i = 2 [c n Cii + c ^ C;t] 



and the random coupling Ju given by 



Jtj = jjlfoj + €flli\, 



where £,• and 77, are random variables which follow a bimodal distribution 
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In the present work the partition function is obtained within the Lagrangian path integral formalism. In that case, 
the spin operators are represented as bilinear combinations of Grassmann fields (tp, if*). Particularly, the 25 -FvH 
model admits only magnetic states. This restriction for the corresponding partition function is obtained by using the 
Kronecker 6 function {6{h^ + «,| - 1) = ■£- L dxt e' Xiinn+na ~^) J3,|29J|. This procedure allows to write the partition 
function for 25 and 45 situations in a compact form as 



Z{y) = e^ f[ 2^ X dXj6 ' y> ZefAy] 



(5) 



and 

Zef/ty) = J DQp'vie^^ 

where s (s = 2 or 4) is the number of states per site, 
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H is chemical potential, /3 - l/T and yj = fifi or yj = ixj for the 45 -FvH and 25 -FvH models, respectively. 
The action A{y] given in Eq. (O can be build using the hamiltonian (Q]i. Thus, we have: 
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The coupling J/j given in Eq. (0 allows to rewrite the random part of the action as: 
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Similarly, the ferromagnetic part of the action is given as 



(15) 



The Hubbard-Stratonovich transformation is used to linearize the first three terms and the first one from the right 
side of Eqs. (TBI) and ( TT5b . respectively. The last terms in the right side of Eqs. (IX4b -(fT3T> vanish in the thermodynamic 
limit. Therefore, Z{y} defined in Eqs. <f5T> - (TT~2T> is obtained as: 
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In Eq. ( fTST l. it has been used the self-averaging property jj Tjjf(£j> Vj)- The average << ... >>f,,, can be achieved by 
using the probability distribution P{%, iff — P{^)P{ifj where P(£) (or P(rj)) is defined in Eq. @. In this work, the SA 
is assumed, in which m - m{r) and q„ = q„(T) for n =1,2 and 3. Therefore, the internal field in Eq.lfTBl is 



hj = V7(//7^i + /£/<72 + Ofa + £7)43) + V7) '«■ 
The integrals over {q,,} and ra in Eq. ( [ToT l can be solved by the saddle point method which gives 
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Figure 1: Phase diagrams T/J versus Jg/J for several values of ' fi/J and T/J = 0. The solid lines correspond to the second order transition while 
the dashed lines indicate first order transition. The inset shows a comparison between the phase diagrams for the 25 -fvH model (dotted lines) and 
the 4S-fvH model with fi = 0. 

The Fourier transform can be used in Eq. ( fTTI i and the functional integrals over the Grassmann fields (cp, if*) performed 
as well as the subsequent sum over Matsubara's frequencies following close Ref. Il5ll to give: 



MQuQi, M) = 2e v [cosh(j) + cosh( ^8A)] 



(23) 



with A = [7(77 + 0(fij + Qi) - JtjQ, - J{Q 2 + J M) 2 + T 2 . 

For the 2S -FvH model the integral over x in Eq. ( [To*] ) is calculated which leads to the same thermodynamics found 
in Ref. 12811 for the quantum vH model. While, in the case of 45 -FvH model, Z{jj) can also be obtained from Eq.dToTi 
which gives for the grand canonical potential Q., the following expression: 



/3Q.=pjQ 2 + ?j-M 2 -/3n 



-«ln 2[cosh(/3//) + cosh(/3 Va)]» 



(24) 



n( 



where it is assumed Qi - Qi - Q- 



2. Results 

In order to investigate the interplay between charge and spin fluctuations in the 45 -FvH model in the presence of 
a transverse magnetic field T, a set of phase diagrams can be built from the numerical solution of the saddle point SG 
order parameter Q and the magnetization M which is obtained from Eq. ( l24l >. It should be remarked that temperature, 
/o (the strength of the ferromagnetic component of Jtj), \x and Y are given in units of J (the strength of the random 
component of Jij). For numerical results, J — 2 is used. In Figs. Q]|4] phase diagrams temperature versus Jq for several 
values of /j and F are displayed. 
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Figure 2: Phase diagrams T/J versus Jg/J for several values of ' fi/J and T/J = 0. The solid lines correspond to the second order transition while 
the dashed lines indicate first order transition. The inset shows a comparison between the phase diagrams for the 25 -fvH model (dotted lines) and 
the 4S-fvH model with /j = 0. 



Fig- E a ) shows the situation corresponding to the half-filling occupation (jj = 0) and r = 0. For that case, the 
phase diagram of the classical vH model is basically recovered. For small Jo, when the temperature is decreased, 
there is first a transition from the PM phase to the SG one which has M — and Q + 0. Then, there is a second 
transition from the SG phase to the mixed (FM+SG) one in which M#0 and Q + 0. For larger J , SG and FM+SG 
phases are replaced by the ferromagnetic (FM) one with M + and Q — 0. In the inset of Fig. []I a X me 45 -FvH 
model phase diagram with /i — is compared to the 25 -FvH one. These two models present qualitatively the same 
results. However, the phase boundaries of the AS -FvH model at the half-filling appear at lower temperatures than the 
25 -FvH model. That is a consequence of the non-magnetic states allowed in the 45 -FvH model. Figs. [lJb)-|TJcl) show 
how the increase of /i can affect phase boundaries even when F = 0. When the average occupation n moves from the 
half-filling (/v = 0) to doubly occupied sites (ji - 1), the transition temperature Tj between paramagnetism and SG 
phase as well as that one (here called T g ) between SG and FM+SG phases are depressed. Ultimately, this depressing 
of Tf and T g as /i increases leads to the complete suppression of SG and FM+SG. Similarly to the FISG model (see 
discussion in Ref. 12111 ). the increase of /i also changes the nature of the phase transitions which appear in Fig. Q] 
That can be seen, for instance, for T g in Fig. |TJb) when /i - 0.25. Furthermore, for// = 1 (see Fig. Q2d)), the nature 
of the Curie temperature T c is also changed. It appears a tricritical point which can be located (see Appendix) at 
jtnc _ 0.759 with J' nc - 2.278. One last consequence due to the increase of// is that larger values of /n are necessary 
to obtain a FM solution. 

In Fig[2j r = 0.5 while /i assumes the same values already used in Fig. Q] Therefore, the results can be compared 
directly to Fig. Q] which allows to understand the combined effects of F and /i. The phase diagram for /u - is 
shown in Fig. [2ja). In this case, the effect of T is to depress both Tf and T g . However, the subsequent increase of 
fi showed in Figs. |2|b)-|2|d) leads to a distinct situation as compared with FigQ] For instance, the FM+SG phase 
is completely suppressed. In contrast, the SG phase is still preserved. Other important difference in comparison to 
the case displayed in Fig. Q]is concerned with the nature of phase boundaries. For instance, the line transition T g 
remains always a second order one. That is not the case for Tf which becomes a first order one as displayed in Fig. 
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Figure 3: Grand canonical potential fi versus T/J for /x/J = 1.0, T/J = 0.5 and several values of Jg. 

[2jd) where // = 1 . In fact, for this value of chemical potential there is a completely new scenario. The existence of 
the tricritical point in the Curie temperature is preserved being located at T' c nc - 0.729. Nevertheless, there is now 
an additional triple point Triple ~ 0.376 with the coexistence of SG, FM and PM phases. Fig. |3]exhibits the grand 
canonical potential fi as a function of temperature for the SG, PM and FE phases near the triple point. The first order 
transitions are located by comparing fi of these phases, where the stable solution is always that one with minimum fi. 
Particularly, the triple point is obtained when fi of the SG, PM and FE phases assume the same values (see Fig. [3jc)) 

In Fig. |4]phase diagrams with F = 1.0 are displayed. In that case, the FM+SG is completely suppressed for// = 
(see Fig. HJa)). The remaining phase diagrams displayed in Figs. |4|b) -|4|d) present only SG, FM and PM solutions 
for the order parameters. When /u = 0.25 and 0.5, the second order line transition Tf is even more diminished when 
compared to the respective previous situations given in Figs. QJ2] For jj - 1, the complex scenario found in Fig. |2jd) 
is repeated in Fig. |3d) with the presence of tricritical and triple points. However, interestingly, Tf, which is now a 
first order line transition, has a small increase as compared with Fig[2|d). 

The effects of the F increasing on the 25'-FvH model can be observed in the inset of Figs. [2|a) and Ufa), which 
exhibit phase diagrams of the 25-FvH model compared to the 45 one with /i = for F = 0.5 and 1.0, respectively. 
As in the 4-S -FvH model, the magnetic phase is also suppressed by F in the 2S -FvH model. In addition, the phase 
boundaries converge to the same values when temperature decreases. 



3. Conclusion 

The present work has studied the fermionic version of the classical van Hemmen model 1U2I1 for SG in the presence 
of a transverse magnetic field F. The goal is to verify how the interplay between charge and spin fluctuations can affect 
the magnetic ordering in this particular model. It should be mentioned that the knowledge of such interplay can also be 
quite useful in further applications for studying the disordered magnetism when other physical processes are present 
(see, for instance, Ref. 111310 . 
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Figure 4: Phase diagrams T/J versus Jo /J for several values of fi/J and T/J = 1.0. Dashed lines indicate first order transition. The inset shows 
phase diagrams for the 25 -fvH model (dotted lines) and the 45 -fvH model with fi = and T = 1 . 



The mentioned fermionic version has been formulated in two versions. In the first one, the so called 25 -FvH, 
the two nonmagnetic eigenstates of the operator S ? (defined in Eq. ( fT2l have been eliminated. Therefore, the corre- 
sponding thermodynamics obtained is exactly the same one obtained by a previous quantum formulation of the van 
Hemmen model given in Ref. 112811 . On the contrary, when the four original eigenstates are preserved, the solution 
for both saddle point order parameters, the SG (Q) and the magnetization (M) one, indicates that not only the phase 
diagram but also the nature of the phase transitions are deeply affected when the chemical potential fi and/or T are 
varied. For instance, the mixed phase FM+SG can be suppressed as well as there is onset of tricritical and triple points 
by increasing ju and/or T. In fact, the role of T as responsible for the suppression of the FM+SG phase has already 
been found in Ref. [28]. Remarkably, in the present problem, the same result can be obtained by varying only yu, or 
in other words, only by changing the average occupation of sites n. Therefore, one can say that the thermodynamics 
derived from van Hemmen model is either sensitive to a spin flipping mechanism provided by F as well a dilution 
one provided by fi. Surely, when both mechanisms are combined, the complexity of the phase diagram is even more 
intense as can be illustrated by a first order phase transition between the PM and SG (or FE) phases with the presence 
of tricritical and triple points. 

To conclude, we have shown in this work that mechanisms such as dilution and quantum spin flipping given by fi 
and T, respectively, are sources of quite non-trivial effects in the fermionic van Hemmen model. More precisely, those 
effects indicate that the existence of magnetic phases and the nature of its boundaries in that model are extremely 
sensitive to the redistribution of charge provided by fi and F. It should be emphasized that thermodynamics in the 
present work could be derived without the use of sophisticated mathematical techniques such as the replica method. 
It is well known that in the case of SG first order phase transitions, the use of this method implies in considerable 
difficulties as, for example, location of first order phases boundary 11221 l23l |24J, |25[ |30i |3lJ . 
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Appendix A. 

The Landau expansion of the grand canonical potential (Eq. (l24l i) in powers of M and Q is given by 

, B 2 J 4 A , , /3 2 J 4 A , 

/3n = n Q +A 2 Q 2 + ^-^Q 4 + B 2 M 2 + '-^ r M 4 (A.l) 

where 

1 A() I 1 A 

1 + coshySrcoshySyU sinhjgr 

Kfi pi A() 

and 

A'o = coshySyL/ + cosh/JT. (A.4) 

The second order phase transition from the paramagnetic phase to the SG one (to the FE phase) occurs when A 2 - 
(B 2 - 0) and A > 0. The tricritical point is located when B 2 - with A = 0. 
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